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Systematic ab initio LDA calculations were performed for all the typical representatives of recently discov- 
ered class of iron based high-temperature superconductors: REOFe(As,P) (RE=La,Ce,Nd,Sm,Tb), Ba2Fe2As, 
(Sr,Ca)FFeAs, Sr4Sc2 06Fe2P2, LiFeAs and Fe(Se,Te). Non-monotonic behavior of total density of states at the 
Fermi level is observed as a function of anion height relative to Fe layer with maximum at about Aza ~1.37A, 
attributed to changing Fe - As (P,Se,Te) hybridization. This leads to a similar dependence of superconducting 
transition temperature Tc as observed in the experiments. The fit of this dependence to elementary BCS 
theory produces semiquantitative agreement with experimental data for Tc for the whole class of iron based 
superconductors. The similar fit to Allen - Dynes formula underestimates Tc in the vicinity of the maxi- 
mum, signifying the possible importance of non - phonon pairing in this region. These results unambiguously 
demonstrate that the main effect of Tc variation between different types of iron based superconductors is due 
to the corresponding variation of the density of states at the Fermi level. 
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Recent discovery of the new class of iron based 
high-temperature superconductors [T] has ignited al- 
^ most unprecedented stream of experimental and the- 
" 'oretical studies (for the review of an early work see 
[21 [3] ) . Despite the immense progress in understanding 
,of these systems, the nature (mechanism) of supercon- 
ducting pairing, as well as the reasons for high values 
of superconducting temperature Tc, are still under de- 
bate. In a recent preprint [J] Mizuguhci et al. have es- 
tablished an interesting anion height dependence of Tc 
for the typical representatives of Fe - based supercon- 
ductors, demonstrating almost symmetric curve with a 
.peak around Aza ^ 1.37A. Below we present an expla- 
nation of this dependence, attributing it to the effect of 
the appropriate variation of the total density of states 
'at the Fermi level within standard BCS-like approach. 

The main representatives of the class of iron (pnic- 
'tides, chalcogenides) based superconductors known at 
the moment are: 
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1. Doped RElll (RE=La,Ce,Pr,Nd,Sm,Tb,Dy) 
with Tc about 25-55 K, with most typi- 
cal representatives such as LaOi_a;Fa;FeAs 
[DElElilllllOlIIIlIiaand LaOi_:,F^FeP [l3] 
with much lower Tc=6.6 K. 

2. Doped A122 (A=Ba,Sr), such as Bai_a;Ka;Fe2As2 
[Himiiniin] and Tc about 38 K. 



3. 



4. 



5. 
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Lii_2;FeAs with Tc=18 K [T8l[T9]. 
(Sr,Ca,Eu)FFeAs [201 HI] with Tc=36 K 
Sr4(Sc,V)206Fe2(P,As)2 with Tc=17 K [21]. 
FeSe^, FeSei_^Te^ with T^ up to 14 K [IS]. 
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There is now a plenty of papers on LDA (local den- 
sity approximation) calculation of the band structure 
of Lain [301 [311 131], LaOFeP [33], RElll series [26], 
BaFe2As2 [HI [Ml ISS], LiFeAs [351 ISS], (Sr,Ca)FFeAs 
[inil^, Sr42622 [3^ and Fe(S,Se,Te) [31]. Below we 
present some of the results of our continued work on 
LDA electronic structure, along the lines of Refs. [IHl 
\T7[ [551 HH], covering all typical representatives of the 
whole class of iron based superconductors and taking 
into account some new structural data. 

Iron based high-Tc superconductors in general have 
tetragonal structure with the space group PA/nmm 
(REllll, LiFeAs, Sr42622, Fe(Te,Se), SrFFeAs) and 
the space group lA/mmm (Bal22). For the PA/nmm 
systems Fe ions occupy positions (2b)(0.75, 0.25, 0.5), 
and anion ions A=P,As,Se,Te - (2c) (0.25,0. 25, Za) and 
for Bal22 Fe(4d) (0.5,0,0.25), As(4e)(0,0,Za). Corre- 
sponding experimental lattice parameters and atomic 
coordinates (used in this work) are collected in Table 1. 

Physically important electronic bands (those which 
cross the Fermi level) are formed by antibonding Fe(3d)- 
A(p) states of FeA4 tetrahedron layer. Here A denotes 
different types of anions: P,As,Se,Te. To calculate elec- 
tronic structure of compounds listed in Table 1 lin- 



1 



2 



E. Z. Kuchinskii, I. A. Nekrasov, M. V. Sadovskii 



earized mufiin-tin orbitals method (LMTO) [JD] with 
defauh settings was employed (except for Relll sys- 
tems, where Re-4f states were taken as a pseudocore 
states). Obtained resuhs are in good agreement with 
other LDA calculations by other authors. 

Motivated by the results of Ref. [4] we present here 
our LDA calculated total density of states lSl{Ep) as a 
function of anion height Aza with respect to Fe layer. 
Corresponding dependence is plotted in Fig. 1 with cir- 
cles. We can see that N(£'i?) has an interesting behavior 
with clear maximum at about Azq ^1.37A (see also Ta- 
ble 2). Such nonmonotonous behavior can be explained 
by hybridization effects. Namely, as a governing struc- 
tural parameter characterizing hybridization strength 
one can chose a-Fe-a angle - an angle between anions 
(a) and Fe within the same tetrahedron. The value of 
the angle corresponding to the strongest hybridization 
is 109.45°, i.e. for an ideal anion tetrahedron with Fe in 
the very center of it. Other crystal structure parame- 
ters which might be marked as importand here such as 
Fe-Fe, Fe-a or a-a distances are not changed very much 
from system to system and do not have any transparent 
dependence of Aza- The values of these distances are 
about following 2.8A, 2.4A and 3.85A with slight low- 
ering for LaOFeP, LiFeAs and Fe(Te,Se) compounds. 

From Table 1 one can see that compounds with high- 
est lSi{Ep) values have the a-Fe-a angle very close to this 
value. Decrease or increase of this angle leads to N{Ep) 
drop from this maximum value. This comes from par- 
tial DOS behavior. The strongest hybridization corre- 
sponds to the strongest bonding- antibonding splitting. 
Since antibonding band DOS grows monotonically with 
binding energy [26l [27l [28l |29] stronger hybridization 
will lead to higher values of ^{Ef)- With lowering of 
hybridization bonding-antibonding splitting goes down 
together with N{Ef). 

The Aza dependence of N{Ep) inevitably leads to 
the corresponding dependence of superconducting crit- 
ical temperature Tc. To estimate this we, first of all, 
use the elementary BCS expression: Tc = 1.14a;De~^/'^, 
where lud is the characteristic frequency of collective ex- 
citations involved in pairing interaction (phonons, spin 
fluctuations etc.), and A — gN{Ep)/2 is the dimen- 
sionless pairing interaction constant {g is the appropri- 
ate dimensional coupling constant). In the following we 
take aj_D=350 K in rough accord with neutron scatter- 
ing experiments on phonon density of states for Lai 11 
[5T] and Bal22[42 systems. We fix g to fit the exper- 
imental value of Tc for Bal22 system since this system 
possesses probably most stable value of Tc (about 38 K) 
with respect to the way of sample preparation and dop- 
ing. Thus we obtain the value of dimensionless coupling 
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Fig. 1. LDA calculated total DOS values N{Ef) (cir- 
cles, right scale) and superconducting transition tem- 
peratures Tc (left scale) obtained from simple BCS 
(stars) and Allen-Dynes expressions (squares), as well 
as experimental Tc values (triangles) versus anion 
height Aza over Fe layer for a number of iron based 
high temperature superconductors. 

constant A=0.43. Then just fixing the value of g as for 
Bal22 we obtain Tc values for all other systems, tak- 
ing into account the appropriate change of the density 
of states. Rather surprisingly we observe almost quan- 
titative agreement with experimental data on Tc (see 
triangles on Fig. 1 and Table 2). Note that we can even 
obtain the right order of Tc values for 1111 systems with 
different rare-earth elements as due to rather small dif- 
ference of corresponding densities of states, which were 
not obtained in our previous work |26j . where we just 
fixed Aza to the only known at that time experimen- 
tal value for LaOFeAs. However, the calculated value 
of Tc for LaOFeAs system is still rather higher than 
most typical experimental value of 26-28K. At the same 
time, the samples of this system obtained via high pres- 
sure synthesis |43l demonstrated much higher values of 
Tc ^ 41 K, which is pretty closer to our calculated val- 
ues. Also the notable deviation of our calculated Tc for 
LiFeAs system may be attributed both to the crude- 
ness of our model (e.g. our use of a single value of 
for all compounds), as well as to probable experimental 
uncertainties of Tc in this system. 

In principle, for the number of systems under con- 
sideration we can obtain even better results if we use 
the multiple band BCS-like approach, along the lines 
of Ref. ^44|. However, to reduce the number of free 
parameters, the multiple band model fit requires addi- 
tional information on the relations between energy gaps 
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on different Fermi surface sfieets (cylinders), wliicli at 
present is only availbale for some of 122 systems. 

It is well known that the elementary BCS-like ex- 
pression for Tc has a tendency to overestimate the role 
of the density of states at the Fermi level. As an al- 
ternative we try the same approach estimating super- 
conducting critical temperature Tc using Allen - Dynes 
interpolation formula (which is probably the best semi- 
analytic expression for Tc in case of electron-phonon 
pairing mechanism, including the strong coupling re- 
gion) [35]: 



ic = exp 

1.20 ^ 



1.04(1 + A) 
A - /i* - 0.62A^* 



(1) 



where 

/i = [1 + (A/Ai) 

/2 — t H , ,9 1 
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Ai = 2.46(1 + 3.8Af*), 

1.82(1 + 6.3Ai*)(w2/w;„), 



and UJ2 =< oj > ' , a;;„ are square root average and 
average logarithm of phonon frequency. Assuming for 
simplicity cj/„ w aj2=35G K, and taking the optimistic 
value of Coulomb pseudopotential /i* = 0, we repeat 
our previous analysis, fixing first A=0.97 for Bal22 and 
then changing only the density of states as obtained in 
our calculations for all other systems. The results for 
Tc obtained in this way from Allen - Dynes expression 
([ij are shown in Fig. 1 by squares (see also Table 2). 

We can see that Allen - Dynes expression produces a 
kind of a lower bound Tc estimate, with obvious deficit 
in Tc values in the vicinity of maximum. This deficit 
may signify the importance of non phonon pairing mech- 
anism to obtain maximal values of Tc in FeAs supercon- 
ductors. However, our main conclusion on important 
correlation of Tc with the values of the density of states 
at the Fermi level remains intact. 

In fact we do not adhere at the moment to any spe- 
cific pairing mechanism. Main objection to electron- 
phonon pairing in iron based superconductors comes 
from microscopic calculations, e.g. those of Ref. [5T| . 
At the same time, there are experiments on isotope ef- 
fect [inmSj which support the importance of this mech- 
anism, though the other isotope experiments |48) pro- 
duce quite opposite picture. 

Our choice of characteristic phonon frequencies in 
the preexponential factor of BCS and Allen - Dynes ex- 
pressions for Tc is used only as a kind of an estimate. 
What is important to us, is the well known fact that the 
dimensionless pairing constant is proportional to the to- 
tal density of states in almost any BCS-like model of 
superconducting pairing, with some additional modifi- 
cations in the case of multiple band models [33]. 



In conclusion, our results show unambiguous correla- 
tion of the values of superconducting Tc and those of the 
total density of electronic states at the Fermi level for 
the whole class of iron based superconductors, thus sup- 
porting the usual BCS-like pairing mechanism in these 
systems. 
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Table 1. Experimental crystallographic data for iron based superconductors. 



System 


AZa, A 


a, A 


c, A 




Za 


^a — Fe — a 


LaOFeP 


1.140 


3.9636 


8.5122 


0.1487 


0.6339 


104.4 




1.200 


4.0160 


15.543 




0.5772 


105.2 




l.OZU 


4.UOOO 




U. 1410 


U.DOlz 


1 nv 
iU r .0 


CeOFeAs 


1.354 


3.9959 


8.6522 


0.1480 


0.6565 


108.4 


SmOFeAs 


1.357 


3.9270 


8.4413 


0.1420 


0.6608 


108.8 


NdOFcAs 


1.367 


3.9476 


8.5446 


0.1440 


0.6600 


110.5 


TbOFeAs 


1.373 


3.8530 


8.2990 


0.1447 


0.6654 


109.7 


SrFFeAs 


1.369 


4.0110 


8.9650 


0.1598 


0.6527 


108.6 


BaFc2As2 


1.371 


3.9090 


13.2122 




0.3538 


109.3 


CaFFeAs 


1.417 


3.8780 


8.5920 


0.1505 


0.6649 


110.4 


LiFeAs 


1.505 


3.7914 


6.3642 


0.8459 


0.2635 


112.7 


Fe(Se,Te) 


1.630 


3.8215 


6.2695 




0.2599 


111.5 



Table 2. LDA total DOS N(Ef), calculated and experimental Tc for iron based superconductors. 



System 




N(Ei,), 
states/cell/eV 


TBCS^ K 


T^D, K 


Tf"P, K 


LaOFeP 


1.140 


2.06 


3.2 


12 


6.6 


Sr4Sc206Fe2P2 


1.200 


3.24 


19 


27 


17 


LaOFeAs 


1.320 


4.13 


36 


37 


28 


CeOFeAs 


1.354 


4.96 


54 


43 


41 


SmOFeAs 


1.357 


4.66 


48 


37 


53 


NdOFeAs 


1.367 


4.78 


50 


44 


54 


TbOFeAs 


1.373 


4.85 


52 


45 


53 


SrFFeAs 


1.369 


4.26 


38 


39 


36 


BaFe2As2 


1.371 


4.22 


38 


38 


38 


CaFFeAs 


1.417 


4.04 


34 


36 


36 


LiFcAs 


1.505 


3.86 


31 


34 


18 


Fe(Se,Te) 


1.630 


2.02 


3 


11 


14 



